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Abstract: 

A theory of the B Krj' decay is proposed. It is based on the Cabbibo favored h —>■ ccs 
process followed by a direct materialization of the cc pair into the rj'. This mechanism 
works due to a non-valence Zweig rule violating c-quark component of the rj', which is 
unique to its very special nature. This non-perturbative "intrinsic charm" content of the 
Tj' is evaluated using the Operator Product Expansion and QCD low energy theorems. 
Our results are consistent with an unexpectedly large Br{B Krj') ~ 7.8 ■ 10~^ recently 
announced by CLEO. 



1 Introduction 



This paper suggests a theory of the B —>■ Kt]' decay which may shed a new hght on prop- 
erties of the 7]' meson. Our study is motivated by recent resuhs of the CLEO collaboration 
which has announced an unexpectedly large branching ratio 

Br{B Kt]') = {7.8tli ± 1.0) ■ 10"^ (1) 

Little thought is needed to realize that this number is in severe contradiction with a 
standard view of the process at the quark level as a decay of the 6-quark into the light 
quarks which could be suggested as soon as the rj' is usually considered to be a SU(3) 
singlet meson made of the u—, d— and s— quarks (see Sect. 2 for more detail). This result 
may not seem too surprising if one remembers the well known fact that the quark content 
of the 7]' is undistinguishable from the gluon one due to the axial anomaly. One the 
other hand, in the weak decay the 6-quark proceeds more strongly to the ccs system due 
to the Cabbibo enhancement of the latter in comparison to the uus state. Since a pair 
of c-quarks can easily convert to gluons, one can suggest the following scenario of the 
B Krj' decay. The 6-quark proceeds into the s-quark and c-quark pair, while the latter 
directly materializes into the 1]' via a non-valence "intrinsic charm" c-quark component 
of the r]' which exists due to virtual cc ^ gluons transitions. An immediate objection 
to this proposal which can come to one's mind is that this process is expected to bring 
a very small contribution to the decay width as soon as it obviously violates the Zweig 
rule. We will argue that though the scenario we suggest is indeed Zweig rule-violating, it 
nevertheless can explain the data. The reason is that we actually deal here with a situation 
where the Zweig rule itself is badly broken down. As will be discussed in detail in Sect. 4, 
both regularities and sources of breaking down the Zweig rule are nowadays well classified 
and studied. In particular, it is 100 % violated for pseudoscalar mesons including the i]'. 
In effect, we find that an extent to which the Zweig rule is broken down in the problem 
at hand suffices to reconcile the theory with the data (|l]). The uniqueness of rj' is in both 
a possibility to evaluate this effect and its very large magnitude. The decay B Krj' 
serves as a probe of the "intrinsic charm" content of the rj' . On the contrary 
to what it may sound, the mechanism of violating the Zweig rule in the r]' is purely 
non-perturbative. To be honest, we have to note that an accuracy of our result is rather 
low, of the order of factor two in the amplitude. It is important, however, that a main 
source of uncertainty in our approach is well localized and related to a poor knowledge of 
a particular vacuum condensate. Therefore the theoretical precision can be considerably 
improved in the future. We should stress that in contrast to a recent proposal on 
importance of the axial anomaly in the closely related inclusive B — > rj'Xg decay, it does 
not play any role in our mechanism. On the contrary, the anomaly is exactly cancelled in 
the Operator Product Expansion (OPE) in powers of for a c-quark bilinear operator 
(see Eqs.( P^jr7| ) below), which is a starting point of our approach to the problem. In a 
sense, we therefore deal with a "post-anomalous" effect which, of course, is suppressed by 
the parameter l/m^. However, in the real world nic ^ 1.25 GeV is not far from a hadronic 
~ 1 GeV scale and, as will be shown below, an effect of the charmed loop is very large 
numerically. On the other hand, methods applied in our study are in close parallelism 
with those developed earlier in a study of the famous U{1) problem (whose key ingredient 
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is just the axial anomaly), and will be explained in the course of our presentation. 

Our strategy consists of a few steps. We start in Sect. 2 with the standard approach to 
the B Krj' decay and demonstrate that its prediction is about two orders of magnitude 
smaller than the experimental number (|l|). We then propose in Sect. 3 an alternative 
gluon mechanism and explain our method for calculation of a crucial quantity of our 
consideration which is the matrix element (0| 07^750! r/'). Using the data ([l|) as an input, 
we calculate an "experimental" value of this matrix element. To calculate the same 
quantity theoretically, we first reduce it by using the OPE to the matrix element of a 
pseudo-scalar three gluon operator {Q\GGG\rj') . The latter object is further related to a 
particular correlation function of gluon currents extending ideas originally suggested by 
Witten 1^ and Veneziano |Q in their approach to the U{1) problem. This correlation 
function is next calculated in Sect. 4 in terms of a vacuum expectation value of the three 
gluon operator (g^G^) by using QCD low energy theorems. We also discuss there physics 
responsible for breaking down the Zweig rule. In Sect. 5 we estimate the latter vacuum 
condensate and finally obtain a theoretical prediction for the matrix element of interest 
(0|c7^75c|?7'). We compare this number with an "experimental" value found from the 
observed data and find a satisfactory agreement between them. This demonstrates 
that the gluon mechanism indeed explains the data with a reservation for uncertainty of 
our results. A final Sect. 6 contains our conclusions. 



2 The standard approach to B ^ Kr]' 

In this section we estimate a width of the B Kr]' decay assuming that the rj' meson 
is made exclusively of light quarks. In this case the relevant terms in the effective non- 
leptonic Hamiltonian are 
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Here Ot are defined as (we use the notations = 7^(1 — 75), i?^ = 7^(1 + 75) ) 
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where j are the color indices and q is any of the the m, (i, s, c quarks. Gfj,u = G'^^j^t'^ and 
are the gluon and photon field strength tensors. Oi and O2 are the tree level opera- 
tors, while Os-e and O-j-iq are the gluon and electroweak penguin operators, respectively. 
Oil, 12 are the magnetic penguins. The Wilson coefficients q = Q(/i) depend on the renor- 
malization scale /x and to the next-to-leading order ^ ( for asijnz) = 0.118, aem(^z) = 
1/128, mt = 176 GeV, /x ~ 5 GeV ) are given by the set g |, 
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-0.3125 



C2 = 1.1502 , C3 = 0.0174 , C4 



-0.0373 
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C5 = 0.0104 , C6 = -0.0459 , c^ = 1.398 ■ 10"^ , (4) 
cg = 3.919 ■ 10"^ , cg = -0.0103 , cio = 1.987 • 10"^ 
cii = -0.299 , ci2 = -0.634 

Introducing the transition form factor 

{Kiq)\-s^^b\Bip + q)) = 2g^/+(mJ,) + P;. {UH') + /-K')) (5) 

and the rj' residue (the chiral hmit nig = is imphed) 

W\(lWq\0) = -^^/vPm , -^fv = (0-5 - 0.8) ■ ^ ~ 0.04 GeV (6) 

and neglecting for the moment mrji,mK in comparison to m^, the magnetic penguins 
Oil, Oi2 and 0(1/ rub, 1/N) terms in the factorized matrix elements of penguin operators, 
we obtain the following estimate for the amplitude of interest (here N stands for the 
number of colors) 
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where, in particular, we have omitted left-right penguin contributions, which are sup- 
pressed by 1/mfe. In our opinion, this procedure is much better than an alternative 
one, where only a subset of l/rrih corrections is retained. Furthermore, it is well known 
that the factorization does not work in non-leptonic B-decays. Effects due to a non- 
factorizability are usually taken into account in a phenomenological manner by the sub- 
stitution C1+C2/N Oi with ai ^ 0.25 obtained by a global fit of the data on non-leptonic 
B-decays p. Using this number in (^, we end up withQ 

Br{B Kr]') ~ 1 ■ 10"^ (8) 

which is by two order of magnitude smaller than the experimental result fl]]. It is easy to 
see that this small value is a consequence of a small residue of the rj' supplemented with 
the Cabbibo suppression of the h u transition. An account of penguin contributions, 
as can be seen from (|^), does not help much. Indeed, neglecting for simplicity the tree 
level 6 — > u transition, we can obtain an estimate for the ratio of the decay width of the 
process of interest to the width of the decay B Kcj): 



T{B^Ki) ^ m-s^,^,s\Q){K\-s^,h\B)\' 



V{B^K<P) m-s^^s\Q){K\s^,h\B)\^ 




2.5-10"^ (9) 



^We disagree with where much larger Br{B Krj') ~ 3 • 10^*^ was proposed. In our opinion, this 
large width came as a result of an incorrect assignment of absorptive parts to matrix elements of penguin 
operators, which, by definition of the OPE, are not there. In particular, it follows (see p. 2187) from 
the formulas given in jioj that this decay width becomes infinite (?) in the chiral SU(3) limit. In fact, 
at the level of penguin contributions the decays B — > Kr]' and B — > Kcj) are just identical and, assuming 
that factorization works reasonably well, it is simply impossible to obtain anything substantially different 
from our estimate (M). 
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where we have used the definition (0|s7^s|O) = e^f^m^ with ~ 240 MeV known 
experimentally from the (f) — > e"'"e~ decay ( stands for the polarization vector of the 
0-meson). As Br{B Kip) ~ 1 ■ 10~^ , we obtain a very small magnitude Br{B 
Kt]') ~ 2.5 • 10"'^ in reasonable agreement with (|]). It is now obvious that corrections 
due to a non-factorizability of penguin operators, magnetic penguins contributions^, as 
well as 0(l/mb, terms which have been omitted in (|^), cannot substantially change 
the estimate (H). We therefore conclude that the image of the rj' meson as the SU(3) 
singlet quark state, made exclusively of the m, rf, s quarks, is not adequate to the problem 
at hand. To avoid possible misunderstanding, we should note that the axial anomaly is 
in fact taken into account in the above mechanism. However, its role there is merely to 
fix the residue of the quark currrent into the rj' . 



3 The gluon mechanism in B ^ Kr}' 

Here we suggest an alternative mechanism for the B Ki]' decay which is based on 
the well known fact that the rj' is a very special meson strongly coupled to the gluons. 
Therefore, the process of interest can be mediated by the 6 — >• c decay followed by a 
conversion of the c-quarks into gluons. This means that the matrix element 

(0|C7^75C|V(P)) = V (10) 

is non-zero due to the c —>■ gluons transitions. Of course, since one deals here with virtual 
c-quarks, this matrix element is suppressed by the factor. On the other hand, the 

c-quark is not very heavy and, taken together with the Cabbibo enhancement of the 
b ^ c transition in comparison to b ^ u, the suggested scenario of the B —>■ Krj' can 
be expected to successfully compete with the standard one described in Sect. 2. Actually, 
this gluon mechanism will be argued to dominate the decay. To get a feeling of how 
large the residue (0) must be in order to explain the data (|1|), we reverse the arguments 
and estimate this quantity "experimentally" under assumption that the proposed gluon 
mechanism exhausts the B Krj' decay. A corresponding number is easy to calculate. 
In the factorization approximation the amplitude takes the form 

M = ^V,,V:Mv\p)\cWc\0){Kiq)\rs^,b\Bip + q)) (11) 

(here ai ^ 0.25 , see Sect. 2). For the B ^ K transition form factor (|]) we use the dipole 
formula 

1 — /m% 



with /+(0) ~ 0.32 , m^, ~ 5 GeV [|T2|, [T3|. Calculating now the decay width, we obtain 



numerically the branching ratio in terms of the residue f^S^ defined in (ITO| 



Ac) 

Br{B Kri) - 3.92 ■ 10"=^ ■ ( | (13) 



^It has been argued that the magnetic penguin operator 0\\ enhances the branching ratio for the 
& — > 50 decay by 20-30 % [Q. We expect a similar (or, anyway, not larger) effect of this operator for the 
B Kt]' decay. 
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which together with the data (0) imphes the "experimental" value (we use the central 
value of the branching ratio (|l|) ) 

4,'^ ~ 140 Mel/ ("exp") (14) 

This number may seem to be too large for the proposed mechanism to work as it is only 
a few times smaller than the analogously normalized residue (0| 07^750] 77c(p)) = ifricPfj. 
with /^^ ~ 400 MeV known from the rjc — > 27 decay. However, as will be argued in the 

(c) 

rest of this paper, the theory is able to produce such a large residue f^. In effect, the 
gluon mechanism completely overplays the standard one by two orders of magnitude in 
the decay width, and reconciles a theoretical prediction for the B Krj' decay with the 
data. 

We now proceed to a theoretical calculation of the rj' residue of the charmed axial 
current ([101) . Making use of the anomaly equation, we obtain from (p!OD 

= J_(o|2mea75C+ ^G^.G^.IV) (15) 
' m^, 47r 

Since the c-quark is heavier then the r/', it cannot contribute the matrix element (|T^) on a 
valence level. It does, however, contribute when propagating in a loop. The c-quark in the 
loop is subject to external gluon fields populating the rj' . A technical tool which allows 
to evaluate the corresponding contribution to the matrix element (|15]) is the Operator 
Product Expansion in inverse powers of the c— quark mass (the heavy quark expansion): 

2m,d75C = -^G,,G,. - -^^^g^ r^^Gl^Gl^Gl^ + ... (16) 

A detailed derivation of Eq.(|l6]) is given in Appendix, while here we restrict ourselves by a 
few comments. The first term in the right hand side of (ITBI) is the usual anomaly term with 



the opposite sign. This sign can be easily understood if one reminds that the anomaly 
term corresponds to a subtraction of the Pauli-Villars regulator from the naive divergence 
2mqqi'-f5q of the axial current ^7^755'. On the other hand, the Pauli-Villars contribution 
is a special case of the heavy quark expansion (piBf ) with the strict limit Mr = 00. The 
cancellation of leading terms ~ GG which do not depend on rric is in agreement with the 
intuitive idea that heavy quarks cannot contribute matrix elements over light hadrons^ 
in the limit mg — > 00. The second term in (|TB]) is the gluon operator of the lowest (after 
the anomaly term) dimension d = 6. We omit one more d = 6 operator {D^G ^u){DaGav) 
which can be related to a four-quark operator using the equation of motion. One can 

^ A well known example is the problem of a light particle mass: in the conformal anomaly equation 

q ^ 

(where is the trace of the momentum- energy tensor and the sum is taken over all quark flavors) the 
heavy quark can only contribute when propagating in a loop, and its contribution cancels a correspond- 
ing contribution of the heavy quark to the /3-function in the second term .Thus, the c-quark does not 
contribute the 77' mass in the limit rric 00. 
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show that a contribution of this operator is suppressed both parametrically 
as) and numerically, 



in and 



Dots in (|16]) stand for higher dimensional operators which we do 
not address here. To this accuracy we therefore obtain 



Ac) 

In' 



3 pabcria 



/^a r^b r^c 



We note that the gluon operator in (p!7D correspond to a normalization point /i^ 



(17) 



We have thus reduced the problem to a calculation of the matrix element of the 
purely gluonic operator (p!?]). Apart from a trivial rescaling of the normalization point, 
this matrix element is essentially defined by a low energy physics on a scale ~ 1 GeV. 
One could therefore think that we did not make any progress at all as matrix elements 
of gluon operators are usually not easy to calculate. A situation with the 1]' is, however, 
exceptional, and the matrix element ( [T7| ) is amenable to a theoretical study. We will 
now show that it can be evaluated making use of the large A^ line of reasoning along 
with the property m^, ~ 1/A^ in close analogy with a way Witten has addressed a very 
similar matrix element (0|GG|?7'). The fact that m^, ~ 1/A^ was established by Witten 
0] a long time ago in connection to the celebrated U{1) problem. Witten's objective was 
to understand within the large N approach how massless quarks are able to bring the 
correlation function of the topological density 



T{p) 



47r 47r 



at zero momentum p — down to zero as required by the chiral anomaly, if quark loops 
are suppressed by a power of A^, and thus apparently do not show up to leading order in 
1/A^. To lowest order in 1/A^ the two-point function ([T8|) is given by sums over one-hadron 
intermediate states 



T{p) 



glueballs ~ 



+ subtractions 



+ 



E 



A^ ' — ' 

' mesons ' ' "n 



(19) 



where M„, a„ and rrin, N ' Cn stand for the masses and residues of nth glueball and 
meson states, respectively. Here a„, = 0{N^), and moreover to lowest order in 1/A^ 
the glueball residues a„ do not depend on whether massless quarks are present in the 
Lagrangian or not. Therefore the first term in ([19|) to leading order in 1/A^ corresponds to 
the two-point function in pure Yang- Mills theory (gluodynamics). The crucial observation 
made by Witten and Veneziano was that while a cancellation of the two terms in 
(p!9|) is not possible at generic 7^ 0, it can happen at = , if there is a single meson 
with m? ~ 1/A^, which would then cancel the whole sum over glueballs together with 
subtraction terms in ([T9|) . Witten and Veneziano have further identified this meson with 
the r]' since the latter is the lightest flavor singlet preudoscalar state in nature. 

The reason we have repeated at length the argument due to Witten and Veneziano is 
that, as is easy to see, it carries practically without a word of alteration over any non- 
diagonal correlation function of the topological density ^Gfj,iyG^i, and arbitrary gluon 
operator with the O quantum numbers. Choosing for such the three gluon operator 
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defining the matrix element (p!7|), we obtain the relation 

-^ J dx{0\T{g-'r^G;,Gl^Gl^{x) ^G,.G,.(0)}|0)yM + 0(1) (20) 

where the subscript YM means that the correlation function in ( pOD refers to pure Yang- 
Mills theory. Its calculation will be addressed in the next section. Here we would like to 
mention that, as follows from (pO|), the residue of interest {0\g^GGG\r]') = 0{N^^/'^). 



4 QCD low energy theorems 



The two-point function (|20|) is a new unknown quantity which has to be evaluated in order 
to estimate the matrix element (|TO|). Analogously to the diagonal correlation function of 



the topological density (|15), it vanishes in the perturbation theory to all orders in a^, 
since the topological density is a total derivative whose matrix elements are all zero at 
the perturbative level. Therefore the correlation function (pOD can only be non-zero due 
to non-perturbative contributions. The goal of this section is to estimate (pO|). 

The idea of how to study correlation functions of gluonic currents with O^'^^^ quantum 
number (like our Eq. (|2^) ) was suggested long ago , and we would like to shortly repeat 
it here. It has been known for a long time that in channels with scalar or pseudoscalar 
quantum numbers leading contributions to correlation functions of composite operators 
are related not to standard vacuum condensates, but rather to the so-called direct in- 
stantons |^. The motivation to introduce this object into the theory was a very strong 
indication that the standard QCD sum rules [0] with standard power corrections due 
to local vacuum condensates are not able to describe the O"*"*^"^ channels. In different 
words, OPE does not reproduce there a scale of phenomena which is exactly known from 
elsewhere. A source of this effect was found in existence of direct instantons. 

A meaning of this object is best explained if one considers first a two-point function 
of (pseudo-) scalar gluon currents at large Euclidean momentum Q"^ . In this case a 
leading non-perturbative contribution is obtained when the momentum Q is transfered 
as a whole to a second vertex by a vacuum field (this is allowed by quantum numbers of 
the current) which therefore must be of small size p ~ l/Q. Such situation corresponds 
to a small coupling regime, in which the quasiclassical approximation becomes accurate. 
The vacuum field is therefore classical; it is the famous BPST instanton |T^. Thus, a 
recipe for calculation of the direct instanton contribution at large is simple: the gluon 
field in the current must be substituted by the instanton. The integral over the instanton 
size is then dominated by small p ~ 1/Q. However, with going down to a resonance 
region Q"^ ^ a few GeV"^, this simple picture breaks down [§] - instantons start to 
interact strongly with each other and large size vacuum fields, and the one instanton (or, 
what is the same, dilute instanton gas) approximation stops making sense. A consistent 
calculation of instanton contribution in this case becomes a complicated problem which 
requires going beyond the dilute instanton gas approximation and taking into account 
instanton interactions e.g. in a form suggested by the instanton liquid vacuum model of 
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Shuryak and Diakonov-Petrov (see and references therein). We shall not proceed 

with this approach which is basically a specific model of the QCD vacuum. Rather, we will 
follow an alternative method which was proposed by Novikov et al. (NSVZ) [^]. It makes 
use of a strong assumption that though a vacuum field transferring a small momentum 
Q resembles only a little the original undeformed BPST instanton, it nevertheless retains 
the (anti-) self-duahty property of the latter 

G^Lv = ^G°^u (21) 

in absence of massless quarks, i.e. in YM theory (where instanton transitions are not 
suppressed by fermion zero modes). This conjecture has been supported by explicit cal- 
culations of next-to-leading (after direct instantons) non-perturbative corrections for two- 
point functions of the scalar and pseudoscalar GG currents at moderate . Up to an 
overall sign, they turn out the same (the direct instanton contributions are identical in 
both channels by definition). Therefore, a further extrapolation of the self-duality selec- 
tion rule (plD to even lower — is expected to be correct at least to a 100 % accuracy. 
In fact, a phenomenologically successful mass formula for the rj' derived in Ref.p] by using 
the selection rule (|2l| ) implies that an actual accuracy of this approximations is of the 
order of 50 %. This mass formula has been obtained^ by relating the residue of the rj' 
to the value of the two-point function of topological density in YM theory at = 0, 
while the latter was evaluated using the selection rule (^Tj) and a low energy theorem (see 
Eg. (^2]) below). As a result, m'^, was found proportional to the gluon condensate {g'^G'^) 
in YM theory. 

Following the same logic, we therefore assume that the self-duality selection rule (|2T|) 



can be also applied to the correlation function of interest (^Of) . If this is the case, the 
value of the latter is fixed by the low energy theorem 0] 

2d 



r n la 

: y rfx(0|T{O(x) -^G\m^)YM = -^{0)ym (22 



Here 0{x) is arbitrary color singlet local operator of canonical dimension d made of 
gluons and h = 11 /3N stands for the first coefficient of the /3-function in pure Yang-Mills 
theory. As a derivation of the fundamental relation (^) is rather simple, for the sake of 
completeness we would like to remind it here. One starts with a redefinition of the gluon 
field 

Gfj,u = QoG^jy (23) 

where Qq is the bare coupling constant of QCD defined at the cut-off scale Mq. Then the 
path integral representation immediately yields the relation 

t J dx{0\T{O{x) G\0)}\0) = -^^^^{O) (24) 



''Somewhat differently from Witten's arguments, this was done without an explicit reference to the 
large N picture, but rather using the fact that the 77' is light on a characteristic mass scale in the 0~ 
channel, and therefore its inclusion is a local ~ effect in the momentum space which must nullify 
the two-point function (p8|) in full QCD. 
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On the other hand, the renormahzabihty and the dimensional transmutation phenomenon 
in a massless theory (either QCD with massless quarks or gluodynamics) ensure that 



to) = const Moexp --^ (25) 

L V mj. 

with the choice b = 11/3A^ — 2/3nj (where rij is a number of flavors) or 6 = 11/3A^, 
respectively. Finally, performing the differentiation yields the low energy theorem (p2D. 
More accurate derivation (which gives the same flnal result) including a regularization of 
ultra-violet divergences in (|2^) can be found in Note that by definition perturbative 
contributions are always subtracted in vacuum condensates like (O). Using now the low 
energy theorem (|2^) for the particular choice O = g^G^, we obtain 

^ l__{g^G^)YM_ .26^ 

This is the main result of this section. Coming back to our original definition (p!0D , 
we see that we have related the residue of the charmed axial current into the rj' with 
apparently completely unrelated quantity which is the value of cubic gluon condensate 
in pure Yang-Mills theory (the matrix element of the topological density is known |p, 
{0\{as/4:7i)GG\ri') ~ 0.04 GeV^). This object will be addressed below, while here we 
would like to end up this section with a discussion of one important conceptual point. It 
is well known that usually non-diagonal transitions between quarks and gluons or between 
quarks of different flavors are suppressed - this is the famous Zweig rule. A most popular 
theoretical explanation of this phenomenon exploits the large N argumentation: in this 
limit all non-diagonal in flavors two-point functions are suppressed by powers of 
relatively to diagonal ones. Thus, at ^ oo any mixing dies off. On the other hand, the 
Zweig rule is strongly violated in the scalar and pseudoscalar O"*"^"-* channels. A well known 
example of this violation is provided by the pseudoscalar meson nonet: while in the vector 
channel the p— and u— mesons are almost pure mass eigenstates of the broken flavor SU(3) 
p ~ {uu — dd) , uj ^ ss and the p — u mixing is small, nothing similar is observed in the 
pseudoscalar nonet. There the t] is predominantly the octet rj ~ {uu + dd — 2ss), which 
means that the mixing is 100 %. The theoretical explanation as to why the Zweig 
rule is violated for quark or gluon currents with the 0^*^"-* quantum numbers is that in 
these channels there are direct instantons which are able to convert quarks into gluons 
and vice versa qq ^ gg at the classical level without any suppression. Literally speaking, 
all factors ensuring a smallness drop out: powers of the coupling constant disappear since 
the instanton fleld is strong, G^^, ~ !/(/, and geometrical loop factors like l/(167r^) do not 
arise because there are no loops. One of the striking examples of such kind is a conversion 
of gluons into photons while naively the amplitude 



in fact it is only 0{a/TT) according to a strict low energy theorem which reads 

^-p^G'l^ikMk,)) = f 
4a, Stt 



m^G'HkMh)) = ^Nnf{Ql)F^jF^J (28) 
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( Here (Qg) is the mean quark electric charge and = /cW^W — fcW^W , i = 1, 2 stands 
for the fields strength of a plane wave.) We would like to note that this process is rather 
similar to our case. The gluons proceed to photons through a loop of the c-quark, but 
the perturbatively expected suppression factor (as/vr)^ does not occur! 

A natural question to ask now is the following: does all this mean that the large N 
picture is strongly violated in the = O^*-"-* channels ? A related question is whether the 
experimental fact that m^i/rn?^ > 1 despite that m^, ~ , = 0{N^) implies such 
a violation. The answer to both questions is no. Moreover, neither the direct instantons, 
nor the low energy theorems are at variance with the large N picture. The point is that, 
as was explained in Ref.|p, N is the dimensionless parameter, and true mass relations 
look rather like m^, = M"^ /N , where is some mass scale. It is usually tacitly assumed 
within the large N reasoning that this mass scale is universal for all hadrons. That this 
is not the case was demonstrated by NSVZ the mass scale is not universal but 
determined by quantum numbers of a channel considered. More concretely, this mass 
is set by a scale at which the asymptotic freedom is violated in the particular channel. 
The latter depends drastically on whether direct instantons are allowed (which is the 
case in the O"'"'-"-' channels) or not. If they are there, an interaction of external current 
with vacuum fields is very strong, and the asymptotic freedom breaks down at very small 
distances, i.e. a characteristic mass scale in this channel is not the typical hadron mass 
~ m^, but rather much higher. It is therefore clear that the second of the questions 
posed above is not properly formulated: the mass of the rj' should be compared not 
with m^, but to a characteristic mass in the 0~ channel which is ~ 15 GeV^ P). This 
fact is the "experimental" evidence that the 1/A^ argumentation is quite accurate for the 
rj' , and 1/N terms, omitted in (]20|), are small in comparison to both explicitly written 
terms. At the same time, the above consideration explains why the effects which we 
discuss are not negligibly small numerically: despite of the fact that the matrix element 
{Q\g'^GGG\r]') ~ A^~^/^, a large dimensional parameter in front of N~^/'^ is able to make 
it large in reality. 

5 How large is (g^G^) in pure gluodynamics ? 

We are now returning to the mainstream of our consideration. Our task has reduced, 
according to Eq. (|26D , to a determination of the cubic gluon condensate {g^G^)YM in pure 
YM theory. Note that this quantity does not have to (and in fact does not) coincide 
with the cubic condensate {g^G^) in the real world. While for the latter there exists a 
semi-phenomenological estimate ||I8| 

{g^G^) = {0.0Q- 0.1) GeV^ (29) 

obtained within the QCD sum rules approach, it is of no direct use in Eq. (|2^) . Unfor- 
tunately, we are not aware of any method (except, probably, the lattice approach) which 
could reliably calculate {g^G^)YM to an accuracy of, say, 20 %. Because of this uncer- 
tainty, we are unable to get a theoretical prediction for Br{B Krj') with precision 
comparable to the experimental one. What instead will be argued in this section is that 
different estimates of the value {g^G^)YM enable one to claim that the large number (|I|) 
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is within the realms of our current understanding of non-perturbative QCD. We beheve 
that this statement is interesting by itself in view of a failure of the standard approach to 
this problem (see Sect. 2). It will be shown below that {g^G^)YM > {g^G^)QCD (this result 
follows directly from the theory and has a status of theorem), and moreover numerically 

{g'G')YM = (0.4 - 1.4) GeV^ (30) 



The first tool we are going to use is again the low energy theorem (|2^). Let us recall 
how a very similar question on a value of the condensate {g'^G'^)YM was addressed in the 
classical paper First of all, we note that this condensate corresponds to an imaginary 
world in which all quarks are very heavy. This world could be obtained from the real 
one when the masses of light u—, d— and s— quarks are smoothly drawn up to some large 
value TUq. By the decoupling theorem, this mass must not be very large - once it reaches 
the confinement scale /x ^ 200 MeV , the heavy quarks decouple and do not influence any 
more the gluon condensate. One the other hand, the confinement scale /i is not too far 
from the s— quark mass ~ 150 MeV . Therefore, one can expect that a value of {g'^G'^) 
in a world with m„ = md = rus = mg ^ 150 MeV would give a reasonable estimate for 
the YM condensate {g'^G'^)YM- In the linear approximation in mu,md we need to know 
the derivative 

^(^G^) = -zJdx (0|T {^G^(0) qq{x)] |0) (31) 

The latter two-point function is fixed by the low energy theorem (^) to be proportional 
to the quark condensate known from elsewhere. Since the quark condensate is negative, 
it follows that {d/ dmq){g'^G'^) > 0. In this way NSVZ have obtained an estimate 

{^G')yM = {2-3){^G')QCD (32) 

TT TT 

As has been argued in the sign of effect {g^G'^)YM > {g^G^)QCD is in perfect agree- 
ment with the instanton picture. Indeed, raising the quark masses diminishes the chiral 
suppression of instantons and therefore increases (G^). 
Proceeding analogously, we write 

J^yG') = -ijdx {Q\T{g'G\Q) qq{x)] |0) (33) 

The difference from the case of quadratic gluon condensate (|3T|) is that the two-point 
function does not coincide with the low energy theorem as it was in Eq.(^). Therefore, 
its value is not known exactly. Nevertheless, the low energy theorem (^2|) can still be used 
to estimate the correlation function (^3]). To this end, consider the relation (p^ in QCD 
with h = 11/3A^ — 2/3n/ for three different operators 



lJdx{0\T[^G^x)^G\0)}\0) = li^G^) 

z|rfa;(0|T|gg(x)^G2(0)}|0) = ^{qq) (34) 
tJdx{0\T[g''G'{x)^G\0)]\0) = f{g'G') 
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We now assume that these low energy theorems are saturated by some effective glueball 
state cr with m„ ~ 1 GeV. It should be stressed that we do not insist on existence of a 
real narrow glueball resonance with such mass. Actually, introducing such a (fictitious?) 
glueball amounts to an effective description of the physics of 0"^ channel. Analogous 
methods have been used in a similar in spirit problem of a strange content of the nucleon 
||T9| which also deals with the channel. Note that the glueball mass drops out in the 



final result (|37D . Introducing the residues 



(OI^GV) = Ai , {0\qqW) = A2 , (OI/GV) = A3 (35) 



we put Eqs.(piD in the form 



ml b ^ 7r 

A1A2 6 



ml b 
A3A1 12 3 3 



{qq) (36) 



m?. b 



Using these equations and assuming the same scalar glueball dominance in (^Sj), we obtain 

d ;.3^3^_ A3A2_ 36 {qq){g'G-') ^3^^ 



dmg ml b {^G"^) 

As (qq) < 0, the sign of the derivative is fixed: {d/ dmq){g^G'^) > 0. Moreover, while an 
expected accuracy of the estimate ( p7| ) is of the order of 100 %, the sign is entirely model 
independent and in fact is fixed by the positivity of spectral densities in Eqs.(0). By the 
same token as for {d/ dmg) {g^G"^) , this sign agrees with expectations based on the instanton 
picture of the vacuum. Numerically, using the valuesQ {qq) ~ —0.017 GeV^, {{as/'K)G'^) ~ 
0.012 GeV^, in the linear approximation in Eg . (|37|) yields for = 150 MeV 

{g^G')yM ^ AA{g^G^)QCD = (0.26 - 0.44) GeV"" (38) 



where we have used the estimate (p9|). Note that the effect of going from QCD to gluody- 



namics is larger for the cubic condensate (pSj) than for the quadratic one (|3^). This fact 



is a direct consequence of the low energy theorem (|3^), and can also be easily understood 
on dimensional grounds: since {g^G^) ~ {g'^G'^)^^'^, increasing (g^G^) by a factor ~ 2.5 
(see Eq.(p^)) yields raising of {g^G^) by (2.5)^/^ ~ 4 in comparison to its value in QCD. 
Eq. ( pSD constitutes our first estimate for a value of the cubic gluon condensate in pure YM 
theory which is based on the low energy theorems ( |5^ and the semi-phenomenological 
information on its value in real world (|2£ 



'''The number for the quark condensate corresponds to the normalization point 1 GeV^. Here we would 
like to point out that the normalization pont in our Eq.(p^) is ~ m^, while the low energy theorems 
refer, strictly speaking, to much lower fj, ~ 500 MeV. In view of a large numerical uncertainty of our 
results we neglect this perturbative evolution. Still, one has to bear in mind that the large anomalous 
dimension of the operator 7 — —18 [ pO| is working in the same direction: taking it into account is 
only able to enhance our final result (pQ) or the estimate (pq) by a factor of two. 
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Another possible source of information on the value {g^G^)YM is the instanton liquid 
vacuum model (see [p!5| , |16[). Two basic parameters of this model are average instanton 



size pc and inter-instanton separation R. The latter parameter is chosen such that to 
reproduce the phenomenological value of the gluon condensate {g^G"^) using the fact that 
each instanton contributes a fixed amount / dx g^G^ = 327r^ to this quantity. This yields 
the number R ^ 1/200 MeV~^ for the phenomenological value (in QCD) {{as/'n')G'^) 
0.012 GeV^. This number provides the upper limit for the instanton density in QCD as 
it implies that the entire gluon condensate is due to instantons. On the other hand, the 
ratio Pc/ R does not depend on the value of condensate and is fixed dynamically to be 
Pc/ R — 1/3. A value of the cubic condensate calculated in this model was found to be 
essentially larger than the semi-phenomenological number (pQ]): 



= ^ - 0.9 GeV , {g'G') ~ 0.4 GeV (39) 
{g^G^} 5pi 



(the formula for (g^G^) in terms of the instanton radius pc was first established in [|14[). 
We would like to make the following comment in reference to the result ([391) . For our 
purposes it is more suitable to discuss the instanton vacuum picture not in QCD, but 
in pure YM theory. We note that the instanton vacuum is more simple in gluodynamics 
than in QCD because of absence of the chiral suppression of instantons. In this case 
the inter-instanton separation must be chosen to fit the gluon condensate {g'^G'^) in YM 
theory, which is according to (|3^) larger than the corresponding number in QCD. On the 
other hand, for gluodynamics the ratio Pc/R — 1/3 remains the same. Using (|39|) and the 
estimate (p^), we obtain 

{g'G')YM^ 1.7 GeV^ (40) 

which is a few times larger than our first estimate (|38D . We feel that this result provides the 
upper estimate for the quantity of interest, and the true answer for (g^G^) lies somewhere 
in between of the two numbers (|38D and (0). 

Finally, we would like to discuss information on vacuum condensates, which is available 
from lattice simulations. The quadratic gluon condensate in YM theory on the lattice was 
reported to be 

' 0.15 GeV^ , SU{2) ^ 
^ 0.10 GeV^ , 5f/(3) |22| 

We do not feel qualified enough to discuss a precision of these calculations. In particular, it 
is not very clear (at least to us) whether the large scale separation is accurately performed 
to make possible a comparison with the SVZ definition. Still, it seems undoubtful that 
these results point in the same direction as (^): the quadratic gluon condensate in 
YM theory is essentially larger than in QCD. As for the cubic condensate, though a 
corresponding Monte Carlo data does exist, it is rather difficult to extract from it a value 



i^G^) 

TT 



of the non-perturbative part of {g G ) |2^. An estimate of this quantity can still be 



obtained in an indirect way as the lattice simulations of instantons suggest^ that the 
average size of instantons in gluodynamics on the lattice is approximately 1/400 MeV 
which is a bit larger than the value predicted by the instanton liquid model. In this case 

^E.V. Shuryak, private communication. See also p6|. 
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Eq. ( pUD together with the above value of the quadratic condensate yields an estimate (for 
the SU(3) color group) 

{g'G'')YM^l.^GeV^ (41) 

which is numerically close to (^). Comparing finally all three estimates (^),(^Dp and 
(PI), we suggest that while (|38|) presumably gives a lower bound for the number of interest, 



(^) and ( PI) seem to set up an upper limit with a possible short distance enhancement. 
A reasonable compromise yields our final estimate given above by Eq. (pOj) . 

We are now in a position to estimate the principal input in (|TT| ) which is the residue 
/^f'* of the charmed axial current into the r/', and which was the main object of our 
consideration in this paper. Using ( |26| ) and (pOj), we obtain the following answer for this 
parameter: 

flf = (50 - 180) MeV (42) 

Note that literally the "experimental" number (|I^) corresponds to the value of the con- 
densate {g^G^)YM — 1 GeV^ which is about a midpoint of our prediction (|30|) . Given 
the accuracy of our result (0), we thus conclude that the gluon mechanism seems to be 
sufficient to describe the data (||). Unfortunately, we are currently unable to improve our 
estimate (|42|) , where the main source of uncertainty is due to a poor knowledge of the 
cubic condensate in YM theory (|30|). Some ways to do this will be discussed in the next 
section. 



6 Conclusions 

In this paper we proposed a theory of the B Kt]' decay. We showed that at the quark 
level this process proceeds via the b — > ccs weak decay followed by a conversion of the 
c-quark pair directly into the r]' which is possible due to a presence of a non- valence Zweig 
rule violating "intrinsic charm" component of the rj' wave function. We have found that 
a mechanism of breaking down the Zweig rule in our case is of a purely non-perturbative 
origin. We have further evaluated a most important ingredient of the factorized B — > Kt]' 
amplitude, which is the matrix element of the charmed current (0|c7^75c|?7'), using a 
combination of the Operator Product Expansion technique, large N approach and QCD 
low energy theorems. Our results demonstrate that the proposed mechanism is likely to 
exhaust an extremely large branching ratio measured by the CLEO collaboration, with a 
certain reservation for a poor accuracy of our final answer (|l2|). We do not pretend to have 
made a numerically reliable calculation of the matrix element (0|c7^75c|r7'(p)) = if^f^Pfi 
which in our opinion determines the B — >■ Krj' decay width. We have rather presented 
a semi-quantitative (but parametrically well established) picture which demonstrates a 
close relation between this matrix element and properties of the vacuum of QCD (or, 
more accurately, YM theory). Though the obtained result looks unexpectedly large, it 
has a good explanation within QCD, and is related to strong fluctuations in vacuum o+(-) 
channels. The main source of uncertainty in our approach is well localized and related 
to a value of the cubic condensate {g^G^)YM which is currently not known with enough 
precision. Yet, the exact low energy theorems indicate that this matrix element is quite 
large. The instanton vacuum model and lattice calculations also seem to favor a large 
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value of this condensate. In this reference, more refined calculations of this quantity are 
highly welcome. In particular, it would be very interesting if this condensate could be 
reliably extracted from lattice simulations. An alternative way which can be suggested 
to improve the determination of the condensate {g^G^)YM is akin to the idea of the QCD 
sum rule approach . If we had other physical processes which essentially depend on the 



same cubic condensate, it could be then fixed "phenomenologically" once and forever with 
an expected consistency between theoretical predictions for different physical amplitudes. 
We plan to return to these issues elsewhere. 

We would like to emphasize that the conclusion on a large Zweig rule violating c-quark 
component of the rj' certainly goes beyond the particular example of the B Krj' decay 
and may well be important in other physical processes. We repeat that there are two basic 
reasons for a large magnitude of the residue f^f^: (1) the c-quark mass is not too far from 
the hadron scale 1 GeV and (2) the Zweig rule is badly broken down in the 0^*^"-' channels. 
While there is nothing particularly special (except for its numerical effect) about the first 
factor, the second one is specific to the unique nature and quantum numbers of the rj'. 
Therefore, we do not expect that any other than rj' light particle could yield a similar 
contribution to the B decay. In a more general content, there is an increasing evidence for 
importance of non-valence Zweig rule violating components in hadrons. We remind that 
the problem of the strange quarks in the nucleon (the so-called vrA^ a— term) is resolved 
| |TP| (see also |^) within physics which is very similar to that discussed in this paper. 
Furthermore, there are many other examples where "intrinsic" non-valence configurations 
(including, in particular, the "intrinsic charm" hadron components) seem important (see 



e.g. a recent review p5[ and ||2^). All these examples unambiguously demonstrate that 
non- valence components of hadron can be sizeable. In QCD terms such a situation means 
that a corresponding matrix element has a non-perturbative origin without the naively 
expected a^/vr suppression. We have shown in this paper that this experimentally testable 
physics is amenable to a theoretical control. In this respect, the rj' from B ^ rj' decays 
is an excellent laboratory for a study of fundamental properties of strongly interacting 
QCD. 
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Appendix 



The purpose of this Appendix is to derive Eq.(0) in the text. A convenient machinery 
for such class of problems was invented by Schwinger a long time ago . The Schwinger 



technique allows in many instances (e.g. when one is interested in a short distance ex- 
pansion) to operate with a propagator in external field without a specification of the 
field. A result has a form of expansion in powers of the field and its derivatives. The 
relevant object in our problem is the c-quark propagator at short distances ~ 1/mc, and 
the expansion in the external gluon field amounts to a representation of the propagator 
in a form of series in powers of Kqcd/^c which is a sort of OPE. We refer the interested 



reader to a pedagogical technical review EH] for more detail and relevant references 



The Schwinger operator approach is based on a realization of commutation relations 
of the coordinate and momentum {P^ = iD^ , where is the covariant derivative) 
operators 

[Pf^,X^] = ig^^ 

[P^,P.] = igG%T'' (A.l) 

where are the generators of the color group. One introduces in the coordinate space a 
formal complete set of states \x) as the eigenstates of the coordinate operator : 

X^j\x) = x^\x) 
{y\x) = 5{x-y) (A.2) 

dx\x){x\ = 1 



while in this basis the momentum operator P^ acts as the covariant derivative 

- -I- 

dx 



{y\P^\x) = (z^ + gAl{x)TA 6{x - y) (A.3) 



J 

In these notations we have to evaluate the expression 

d75C = (x|Tr |75(/'— m)~^| \x) (A. 4) 

where m is the c-quark mass and Tr denotes the trace over both color and Lorentz indices. 
Using a resolution of unity 

1 = (f +m)"^(f +m) 

and the formula ( a^^^ = i/2(7^7 



we expand ( |A.4| ) in powers of aG = a^^G'^^JT'^ 



ci^^c = m(x|Tr|75^^^5^^^^(aG)-^5^^--^(aG) 

- V rP2 ^ 2^2 ^^G)—^^{aG)—^^{aG) + --]\x) (A.5) 
8 (P^ — m^j"^ P^ — m'^ P"^ — m'^ J 

= D1 + D2 + ... 
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where dots stand for higher power terms in the external field expansion. It is convenient 
to start the calculation with the second term in ( [A.5|) . It explicitly contains the gluon 
field Gf^u to the third power, and therefore to our accuracy one can neglect the non- 
commutativity of the operators Pa and G^^. Thus 

= -C-^-{(75^^^75(-G)^| (A.6) 

To the same accuracy the momentum operator in the denominator of (|A.6| ) is substi- 
tuted by the c-number p^, and the x-integration is performing using the formula 



which yields 



(27r)4 ' 16TTH{n - 2){n ~ 3) 



28 ■ 3 ■ 7r2 ■ m3 

A calculation of the trace over the Lorentz indices gives 



Tr {i^icrGf} = -2^ ■ tr.G^.G'.^G^^ 
where tVc stands for the trace over the color indices. The latter is 



J. rabc/^a /^b /~ic _ (~^(~^ 

and finally we obtain 

^2 = ^ \ i9'GGG (A.8) 

A calculation of the first term in (|A.5|) is more tedious as now we have to take into 
account the non-commutativity of operators in order to evaluate this expression to the 
GGG accuracy. Using the identity 



2 



P2 — 77^2 p2 _ ^2 p2 _ ^ 

we write 

^1 = T"^<^l^^ {^^ (p2 + 75 (p, _'^2)3 [^^ ("^)]p^("^)} '^^ 

= II1 + U2 (A. 10) 
The first term in (|A. 10|) is readily calculated to yield 

^1 ^ ~Qo — 2 G^uGfj_y , (A. 11) 

62 ■ ■ m 
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while the second one needs more care. Let us perform first the trace over the Lorentz 
indices: 

n, = -2^/(a:|Tr| ^^^_^^^^3 [P^G,.]p^G,.| \x) (A.12) 
Calculating the commutator 

[P\ G,,] = Px[Px, G^,] + [Pa, G,,]Px = 2iPx{DxG,,) + D^G,, 

we obtain 

- 2ig^m{x\Tr | _^^2)3 (-P'^A^'-) p2 ^ ^2 ^m^} I^) 

+ Vm(x|rr I _^^,^J P^PApAG,.)]-^^^^G,.| |a:) 

- 2ig^m{x\Tr | _^ ^2)4 (-P'^a^'-)^a^'-} I^) + ■ ■ ■ (A.13) 

One can easily see that the first term in this expression gives rise to higher dimensional 
operators and thus does not contribute to our accuracy. However, the second term does 
contain a needed power of the gluon field since 

[P\ PxiD^G,,)] = 2tPxPpDpDxG,, + ...= z{Pa, P^jD^D^G^, + ... 

which results in 

n2 = 2tg^m{x\Trl^j^^^—^{Px,Pf3}Df,DxG^,G^, - jpT^—^{D^G^,,)G^Jj \x) 

(A.14) 

In the first term to our accuracy one can substitute P^ ^ Pfj.- We further use the Bianki 
identity to evaluate 

{D'^G^u)G^jy = Da{—D^Gai_i — DfiG,^a)G^i, = 2iGafiGfiuG^a (A. 15) 

and eventually obtain 

"'- = - 32.3.l^..„3 «'ggg 

Finally, collecting together (|A.8| ), ( |A.11|) and ( |A.16 ) and multiplying the whole answer 
by 2m, we arrive at 

2mci-f5C = -—GG ]-^g^GGG+... (A.17) 

47r IGn'^m^ 



which completes the proof of Eg. (]T6|) . 
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